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ABSTRACT 

' The discovery of OGLE 2005-BLG-390Lb, the first cool rocky/icy exoplanet, impressively 

demonstrated the sensitivity of the microlensing technique to extra-solar planets below 

1 1 ■ 10 M®. A planet of 1 M® instead of the expected 5 M® for OGLE 2005-BLG-390Lb (with 

! \ an uncertainty factor of two) in the same spot would have provided a detectable deviation 

with an amplitude of ~ 3 per cent and a duration of ~ 12 h. While a standard sampling 
' interval of 1.5 to 2.5 hours for microlensing follow-up observations appears to be insuffi- 

\Q [ cient for characterizing such light curve anomalies and thereby claiming the discovery of the 

I/"") . planets that caused these, an early detection of a deviation could trigger higher-cadence sam- 

CN| ' pling which would have allowed the discovery of an Earth-mass planet in this case. Here, we 

, describe the implementation of an automated anomaly detector, embedded into the eSTAR 

■ system, that profits from immediate feedback provided by the robotic telescopes that form 

£ — . ' the RoboNet-1.0 network. It went into operation for the 2007 microlensing observing season. 

C^> . As part of our discussion about an optimal strategy for planet detection, we shed some new 

ILJ ' light on whether concentrating on highly-magnified events is promising and planets in the 

. £^ \ 'resonant' angular separation equal to the angular Einstein radius are revealed most easily. 

. Given that sub-Neptune mass planets can be considered being common around the host stars 

$H ' probed by microlensing (preferentially M- and K-dwarfs), the higher number of events that 

can be monitored with a network of 2m telescopes and the increased detection efficiency for 
planets below 5 M® arising from an optimized strategy gives a common effort of current mi- 
crolensing campaigns a fair chance to detect an Earth-mass planet (from the ground) ahead 
of the COROT or Kepler missions. The detection limit of gravitational microlensing extends 
even below 0.1 M®, but such planets are not very likely to be detected from current cam- 
paigns. However, these will be within the reach of high-cadence monitoring with a network 
of wide-field telescopes or a space-based telescope. 

Key words: planetary systems - gravitational lensing - methods: observational. 



1 INTRODUCTION 

After iMao & Paczvnskil J 199 lh first pointed out that microlens- 
ing events can be used to infer the presence of extra-solar plan- 
ets or place limits on their abundance, this techniq ue has now 
become established with several claimed detections i Bond et al.l 



l2004IUdalskietalj20 05: Beaulieu et al. 2006; Gould et al. 2006) 
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The discovery of OG LE 2005-BLG-390Lb teeaulieu et ai]|2006t 
iDominik et"afl l2006). estimated to be 5 times more massive than 
Earth, with an uncertainty factor of two, under the lead of the 
PLANET (Probing Lensing Anomalies NETwork)/RoboNet cam- 
paign demonstrated that microlensing not only can detect massive 
gas giants, but also planets that harbour a rocky/icy surface under a 
thin atmosphere. Moreover, it provided the first observational hint 
that cool rocky/icy planets are actually quite common, as previ- 
ously predicted b y simulations ba sed on core-accretion models of 



ously predicted by simulations ba; 
planet formation dlda & Lir]|2005h 
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It was already estimated bv lBennett & Rhid Jl996h that there 
is a non-negligible chance of 1-2 per cent for detecting an Earth- 
mass planet located at about 2 AU from its host star by means of 
observing a few-per-cent deviation in a microlensing light curve. 
However, such a discovery requires photometric measurements on 
a few hundred microlensing events, assuming that a fair fraction of 
the host stars are orbited by such planets. 

A sufficient number of events can only arise from monitoring 
dense fields of stars. With a probability of ~ 10 -6 for a star in the 
Galactic bulge being magnified by more than 34 per cent at any 
given time due to the bending of light caused by the gravitationa l 
field of an intervening foreground star (Kiraga & Paczynski 1994), 
and such a microlensing event lasting of the order of a month, one 
namely needs to monitor 10 7 to 10 8 stars. This was achieved by mi- 
crolensing surveys like OGLE (Optical Gravitational Lensing Ex- 
periment) (Udal ski et al J 19921) . MACHO (MAssive Compact Halo 
Objects) JAlcock et al . 1993), EROS (Experience de la Recherche 
d'Objets Sombres) JAubourg et al.ll 199 3) and MOA (Microlensing 
Observations in Astrophysics) dMuraki et al.ll 1999 ) with a roughly 
daily sampling. Moreover, all these surveys have been equipped 
with real-time alert systems judalski et all 1 1994 lud alski 2003; 
lAlcock et al1ll99dlGlicensteinll200ll ; lBond et alj|200ll) that notify 
the scientific community about ongoing microlensing events. This 
allows to schedule follow-up observations that provide an increased 
photometric accuracy, a denser event sampling, and/or coverage 
during epochs outside the target visibility from the telescope site 
used by the respective survey campaign. 

The PLANET (Probing Lensing Anomalies NETwork) col- 
laboratiorQ established the first telescope network capable of 
round-the-clock nearl y-continuous hig h-precision monitoring of 
microlensing events iAlbrowetal] [l998) with the goal to detect 
gas giant planets and to determine their abundance. For being able 
to detect deviations of 5 per cent, PLANET aims at a 1-2 per cent 
photometric accuracy. With a typical sampling interval of 1.5 to 2.5 
hrs allowing a characterization of planetary anomalies on the basis 
of at least 10-15 data points taken while these last, the required ex- 
posure time then limits the number of events that can be monitored. 
For bright (giant) stars, exposure times of a few minutes are suffi- 
cient, so that PLANET can monitor about 20 events each night or 
75 events per observing season, but this reduces to about 6 events 
each night or 20 events p er season for fainter st ars, for which expo- 
sure times reach 20 min jDominik et aT1l2002l) . In 1999, MACHO 
and OGLE-II together provided about 100 microlensing alerts, out 
of which only 7 were on giant source stars. This severely limited 
PLANET in its planet detection capabilities: rather than 75 events, 
only about 25 could be monitored per season. The OGLE-III up- 
grade, in effect from 2002, had a major impact on the potential 
of microlensing planet searches, paving the way towards the now 
nearly 1000 microlensing events per year provided by the alert sys- 
tems of the OGLB3 and MOaEj surveys. The much larger number 
of events arising from this upgrade allowed OGLE i tself to obtain 
mean i ngful constraints on p lanets of Jupiter mass jTsapras et al.l 
120031 : ISnodgrass et aljpiool) . while OGLE and MOA have even 
demo nstrated that such planets can in fact be detected by their sur- 
veys l lBond et ai]|2004h . However, for studying less massive plan- 
ets, their sampling is insufficient. At the same time, the OGLE-III 
upgrade enabled PLANET to exploit its full theoretical capability, 



and moreover, it gave PLANET a reliable chance to detect planets 
of a few Earth masses provided that these are not rare around the 
stars that cause the microlensing events . The discovery of O GLE 
2005-BLG-390Lb teeaulieu et aUl2006l : iDominik et alj|2006h ex- 
plicitly proved the sensitivity of the PLANET observations to plan- 
ets in that mass range. 

Microlensing events are also regularly monitored by the 
MicroFUN (Microlensing Follow-Up Network) tearrQ. However, 
rather than exploiting a permanent network, MicroFUN concen- 
trates on particularly promising events and activates target-of- 
opportunity observations should such an event be in progress. Be- 
sides lm-class telescopes, their stand-by network includes a larger 
number of small (down to 0.3m diameter) telescopes operated by 
amateur astronomers, which are well suited to observe the peaks of 
events over which the source star makes a bright target. 

Since the PLANET network is restricted in its capabilities 
of monitoring ~25 per cent of the currently alerted events with 
the observational requirements, the planet detection rate could be 
boosted by using larger (2m) telescopes or clusters of lm-class 
telescopes. In fact, such an upgrade is required in order to ob- 
tain a sample that allows a reliable test of models of the forma- 
tion a nd evolution of plan ets around K- and M-dwarfs. RoboNet- 
1 .fjf| ( Burgd orf et al.l2007h marks the prototype of a network of 2m 
robotic telescopes, not only allowing a fast response time, but also 
a flexible scheduling by means o f the multi-agent contract model 
provided by the e STAR projec{f] dAllan. Na vlor & Saunders ! 20061 : 
Allan et al.j|2006f) . eSTAR is a key player in the Heterogeneous 
Telescope Networks (HTN) consortium and involved in the IVOA 
(International Virtual Observatory Alliance) standards process. 

If one aims at the discovery of Earth-mass planets, the stan- 
dard follow-up sampling of 1.5 hrs usually does not produce the 
amount of data required to characterize the corresponding signals, 
and with less frequent sampling one even faces a significant risk 
of missing any hint for a deviation from an ordinary microlensing 
light curve. However, planets of Earth mass and even below can be 
discovered by shortening the sampling interval to ~ 10 min once a 
regularly sampled point is suspected to depart from a model light 
curve that represents a system without planet. In order to properly 
trigger such anomaly alerts, all incoming data need to be checked 
immediately, and prompt action needs to be taken within less than 
~ 15 min. The amount of data and the required response time for 
achieving a good detection efficiency for Earth-mass planets are 
however prohibitive for relying on human inspection. Therefore, 
we here describe the implementation of an automated anomaly de- 
tector that exploits the opportunities of immediate response and 
flexible scheduling of a network of robotic telescopes. A first sim- 
ilar warnin g system, d ubbed EEWS, had been installed by OGLE 
in 2003 (Udalski 2003), which however involves further human in- 
spection and operates with a single telescope. In contrast, our de- 
sign needs to succeed without any human intervention and take 
care of a heterogeneous telescope network. The underlying algo- 
rithm follows previous experience on the assessment of anomalies. 
We explicitly aim at reaching a significant detection efficiency to 
Earth-mass planets with the current survey/follow-up strategy of 
microlensing planet searches. 

This paper is organized as follows. In Sect.[2]we describe the 
modelling of ordinary microlensing events with particular empha- 
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sis on the importance of robust parameter estimates, not confused 
by outliers, in order to properly identify real deviations. While 
Sect.[5]deals with the general strategy for detecting low-mass plan- 
ets by microlensing, we derive a suitable concept for an anomaly 
detector in Sect. [4] The embedding of the SIGNALMEN anomaly 
detector, that went into operation for the 2007 microlensing cam- 
paign, into the eSTAR project is discussed in Sect. [5] before its al- 
gorithm is described in Sect.[6] Sect.|7]then discusses the prospects 
of the SIGNALMEN anomaly detector for discovering planets of 
Earth mass and below. In Sect. [8] we provide a short summary and 
final conclusions. 

The Appendix makes a point on the inability to detect planets 
at the resonant separation in some of the observed events that was 
not discussed earlier. 



2 ORDINARY LIGHT CURVES AND ANOMALIES 

The bending of light due to the gravitational field of a foreground 
'lens' star with mass M at distance Dl causes an observ ed back- 
groun d 'source' star at distance Ds to be magnified by teinsteirj 



A(u) = 



vf- + 2 



iVu 2 +4 



(1) 



if both objects are separated on the sky by the angle u 6e with 6e 
denoting the angular Einstein radius 



AGM 



(2) 



With the assumption that lens and source star move uniformly, 
where n is the absolute value of their relative proper motion, the 
separation angle can be parametrized as 



u(t) 



< + 



(3) 



with uo denotes the closest approach at epoch to, and tE 
is a characteristic event time-scale. 

Each set of observations with a specific telescope and filter 
comprises a data archive s of observed fluxes and their error 
bars a F . [ s ] at epochs t ' s ' . With the source flux and background 

flux F^ depending on the data archive s, one observes symmetric 
light curves 



F^(t) = F^Alu(t)]+F ] 
peaking at to- 



(4) 



Estimates for (to, tE, uq, Fg , F^) can then be obtained by 



minimizing 



m n l 1 

EE 

= l i=l 



F^(t)-Fl h] 



(5) 



While we use the CERN library routine MINUIT for determining 
(toj tE) ^o)> the source and background fluxes Fj, 3 ^ and F^ for any 
choice of (to, tE, ito) simply follow from linear regression as 



F* = 



F, 



<L > 2^7? ~ ~7f 



F B = 



E 



— 7^~ 



7? ~ 2^ ~7^~ 1^ 



7i — 



E 



— 7^- 



■e^-(e 



A(tj) 



(6) 



where the summations run from 1 to n' fc ', <Ji = ap,, and the in- 
dex [k] has been dropped. Any archive s can only be included if it 
contains at least 3 data points. 

The characteristic form of the light curve described by Eq. 10 
is based on the assumption that both source and lens star are sin- 
gle point-like objects that are moving uniformly with respect to 
each other as seen from Earth. Apart from planets orbiting the lens 
star, significant deviations, so-called anomalies, can however also 
be caused by binarity or multiplicity of lens or source, the finite 
angular size of the stars, or the revolution of the Earth (parallax 
effect). 

Since it is our primary goal to detect light curve anomalies, 
it is essential to ensure that our adopted model is reasonably cor- 
rect. However, frequently our data do not allow strong constraints 
to be placed on the model, in particular during early phases of the 
event. It is a well-known fact that OGLE announce a fair fraction 
of their events with the prediction of quite high peak magnification, 
whereas it turns out later that most of thes e peak at muc h lower 
magnifications. As studied in some detail bv lAlbrowl ( l2004h . this is 
related to the fact that \ 2 -minimization is equivalent to obtaining a 
maximum-likelihood estimate of the model parameters if the data 
are assumed to follow a Gaussian distribution, which is biased, i.e. 
its expectation value does not coincide with the true expectation 
value of the considered quantity. Using the statistics of previously 
observed OGLE events, a Bayesian estimate that can be obtained 
by adding an effe ctive penalty f unction to \ 2 comes closer to the 
expectation value dAlbrowl20 04). While the estimated value can be 
tuned by this, one does not fully get around the problem of large 
indeterminacy of the model parameters. 

A further problem arises from the necessity to avoid that our 
model is driven towards data outliers. Otherwise, real anomalies 
would be missed while points matching an ordinary light curve 
would seem deviant. As a consequence, we would face the problem 
of not being able to distinguish between ongoing anomalies and fur- 
ther data requiring an adjustment of model parameters. Therefore, 
we apply a more sophisticated algorithm for estimating the model 
parameters that is rather invulnerable to outliers. 

The model can be made to follow the bulk of the data by 
downweighting points according to their respective residual (e.g. 
Hoagl in. Mosteller & Tukevll 19831) as follows. With the residuals 



(7) 



and the median of their absolute values f ' fe ' for each data archive, 
we give further (bi-square) weight 

2 



[fc] 



1 



Ik] 

r . 




for 
for 



< 



(8) 



to each data point, where we adopt K = 6 for the tuning constant. 
The choice of the weig hts, Eq. l[8), means that data points whose 
absolute residuals exceeds K times their median are ignored. This 
procedure is repeated until the formal x 2 converges. However, we 
need to deal with non-linear models which are prone to several pos- 
sible x 2 minima. In contrast to linear models, it can therefore hap- 
pen that this procedure leads to periodic switching between differ- 
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ent minima, where nevertheless a subsequence converges to each 
of these. In this case, we have to live with the absence of a unique 
minimum and choose that one with the lowest \ 2 - With the formal 
y 2 not being dominated by outliers, we can also reliably adjust the 
relative weight between different data archives k after each itera- 
tion step, so that all (x 2 )' fc V n ' fe ' coincide, preventing the estima- 
tion of model parameters being influenced by the collective over- 
or underestimation of error bars. 



3 DETECTION OF LOW-MASS PLANETS 



It was pointed out bv lMao & Paczynski| dl99 ll) that planets orbit- 
ing the lens star can reveal their existence by causing significant 
deviations to microlensing light curves. They also found that the 
probability to detect a planet becomes resonant if the angular sep- 
aration from its host star is comparable to the angular Einstein 
radius #e, which reflects the fact that the detection of planets is 
aided by the tidal field of their host star. However, as pointed out 
in the Appendix, for a given event, in particular for larger im- 
pact parameters, the detection probability of smaller planets can 
actually drop to zero for angular separations close to #e rather 
than reaching a maximum. In such case, only slightly wider or 
closer separations can be probed. It is a lucky coincidence that the 
gravitational radius of stars and distances within the Milky Way 
combine in such a way that the angular Einstein radius converts 
to a projected separation Dl 9e ~ 2 AU for M — 0.3 M®, 
the typical mass of the lens star s , assu ming Ds ~ 8.5 kpc and 
Z)l ~ 6.5 kpc. iGould & Loebl ( 119920 quantified the prospects 
for detecting planets from microlensing signatures by finding that 
Jupiter-mass planets distributed uniformly within angular separa- 
tions 0.6 (9e ^ d0E ^ 1.6 9e, comprising the so-called leasing 
zone, have a probability of 15 per cent of being detected among 
microlensing events with peak magnifications Ao ^ 1.34, corre- 
sponding to the source entering the Einstein ring (of angular radius 
9e) of the lens star, i.e. no ^ 1. As shown by Griest & Safizadef] 
(1998), this probability increases significantly if one restricts the 
attention to events with larger peak magnifications, where about 
80 per cent is reached for Ao ^ 10. Since the area subtended on 
the sky by angular source positions that correspond to a significant 
deviation decreases towards smaller planet masses, both a shorter 
duration of the planetary signal and a smaller probability to observe 
it result. In contrast, the signal amplitude is only limited by the fi- 
nite angular size of the source, where significant signal reductions 
start arising once it becomes comparable or larger than the size of 
the regio n for which a point sour ce provides a significant deviation. 
However. lBennett & Rhid (l996) estimated that Earth-mass planets 
still have a 1-2 per cent chance of providing a signal in excess of a 
few per cent. 

Planets around the lens star affect the light curve only by 
means of two dimensionless parameters, namely the planet-to-star 
mass ratio q and the separation parameter d, where d 9e is the in- 
stantaneous angular separation of the planets from its host star (i.e. 
the lens star). With typical relative proper motions between lens 
and source stars of /i ~ 15 |i.asd _1 , microlensing events on Galac- 
tic bulge stars are usually observable for about a month or two, 
whereas planetary deviations last between a few hours and a few 
days, depending on the mass of the planet. In contrast to other in- 
direct techniques, microlensing therefore obtains a snapshot mea- 
surement of the planet rather than having to wait for it to complete 
its orbit. This gives microlensing the unique capability of probing 



planets in wide orbits whose periods otherwise easily exceed the 
life-time of a project or its investigator. 

With many events on offer from the OGLE and MOA sur- 
veys and only limited resources available for follow-up observa- 
tions, one needs to make a choice which of these to monitor and 
how frequently to sample each event. With the goal to maximize 
the number of detections of planetary deviations, a prioritization al- 
gorithm that spreads the av ailable observ ing time over the potential 
targets has been devised bv lHornd l !2007h . which forms a central en- 
gine of the RoboNet observing strategy. Any such strategy must be 
based on observables, model para meters arising from the collected 
data, or any other data statistics. As lHornd ( 120071) pointed out, each 
data point carries a detection zone with it, composed of the angular 
positions for which a planet would have caused a detectable devia- 
tion. Unless fi nite- source effects begin diminishing the delegabil- 
ity of planets dHanll2007h . detection zones grow with the current 
magnification. Moreover, the same photometric accuracy can be 
achieved with smaller exposure times for brighter targets. An ef- 
ficient prioritization algorithm therefore needs to be based on both 
the current magnification and brightness along with the time when 
the last observation was carried out, where taking into account the 
latter avoids obtaining redundant information. Such a prioritization 
of events however does not consider how well an observed devia- 
tion allows to constrain its nature of origin and it also assumes that 
the model parameters of the ordinary light curve are known exactly. 

If the effect on the microlensing light curve is dominated by 
a single planet, the lens system can be fairly approximated as a 
binary system consisting of the star and this planet. Gravitational 
lensing by a binary poin t-mass lens has been stu died in great de- 
tail for equal masses by I Schneider & Weifi ( 1986) an d later gen- 
eralized for arbi trary mass ratios by Erdl & Schneider! Jl993l) . On 
the other hand. lChang & Re fsdal ( 1979) have discussed lensing by 
bodies of different mass scales. While their target of interest was 
the brightness variation of individual images of QSOs that are grav- 
itationally lensed by an intervening galaxy, a very similar situation 
arises for planets orbiting a lens star. Similarly to individual stars 
in the galaxy splitting an image due to lensing by the galaxy as 
a whole into 'micro-lensing', a planet can further split one of the 
two images due to lensing by its hos t star if it roughly coincides 
in angular position with that image. iDominikl ( 119991) has further 
investigated the transition towar ds extreme mass rat ios and shown 
how the case described bv lChang & Refsdall h979» . the so-called 
Chang-Refsdal lens, is appr oached . The derived expansions into se- 
ries have later been used by |Bozzal ll999) for discussing the case of 
multiple planets. Binary lenses in general and planetary systems in 
particular create a system of extended caustics, consisting of the an- 
gular positions for which a point-like source star would be infinitely 
magnified. While sufficiently small sources passing the caustics can 
provide quite spectacular signals, planets are more likely to already 
reveal their existence on entering a much larger region surrounding 
these. 

For less massive planets, there are usually two separate regions 
for positions of the source star that lead to detectable planetary sig- 
nals, which are related to two types of caustics. Only if the angular 
separation of the planet from its host star is in a close vicinity to 
the angular Einstein radius 9e, where the corresponding range is 
broader for more massive planets, a single caustic results and these 
regions merge. Otherwise, there are one or two planetary caustics 
which are located around positions for which bending of its light 
due to the gravitational field of the lens star causes the source to 
have image at the position of the planet, and a central caustic which 
can be found near the lens star dGriest & SafizaderJll998l : lDominikl 
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Figure 1. Model light curve of microlensing event OGLE 2005-BLG-390 
along with data taken with the Danish 1.54m at ESO LaSilla (Chile), red, 
the Perth 0.6m (Western Australia), blue, and the Canopus 1.0m (Tas- 
mania), cyan, by PLANET, the Faulkes North 2.0m (Hawaii), green, by 
RoboNet- 1.0, the OGLE 1.3m (Chile), black, and the MOA 0.6m (New 
Zealand), brown, where Am = 2.5 Ig A(t) has been plotted along with 
mi = 2.5 lg A,. The ~ 15 per cent deviation lasting about a day re- 
vealed the existence of a planet with m ~ 5.5 Mgj (uncertain to a factor 
two), while an Earth-mass planet in the same spot would have caused a 3 
per cent deviation lasting about 12 hours (thin line). The time-scale of this 
event is = 11.0 d, while d = 1.610 and q = 7.6 X 10~ 5 . Moreover, 
no = 0.359, to = 31.231 July 2005 UT, and the angle between the vector 
from the planet to its host star and the source trajectory is a = 157.9°, 
where the less centre of mass is to the right hand side. Finally, the source 
star moves by its own radius relative to the lens within t+ = 0.282 d. The 
dotted line refers to a model light curve in the absence of a planet. 

1999). As iBozzaU 19991) demonstrated, the planetary caustics asso- 
ciated with different planets are almost always separated and any 
kind of interferen ce between these is quite unlikely. In contrast, 
lGaudietal.l ( fl998l) pointed out that the central caustic is always af- 
fected by the combined action of all planets. However, it is likely, 
although not guaranteed, that there is a hierarchical order among 
the effects of diffe rent planets, so that a linear super position is a 
fair approximation dRattenburv et alj2 002: Han 2005). 

While the absence of any deviations near the peak of extreme 
highly-magnified ordinary events that are related to the source po- 
tentially approaching the centr al caustic poses strict limits on th e 
abundance of low-mass planets dAbeetal .120 04: Dong et al. 2006), 
their actual discovery from this kind of deviations suffers from 
several complications. While the linear size of the detection re- 
gion around planetary caustics scales with the square root of the 
planet mass, it is proportional to the planet mass itself for the 
central caustic dChang & RefsdaJI 1979): iGriest & Safizadebll 19981 : 
lDominiklll999l : IChung et alj2005l : lHan1(2006t) . Therefore, the finite 
angular size of the source star is more likely to cause a significant 
reduction of the signal amplitude. Moreover, the characterization 
of the nature of origin for such deviations is significantly more dif- 
ficult than for deviations related to planetary caustics. The latter 
provide further information by means of the time elapsed between 
the peak of the background ordinary light curve and the deviation, 
whereas central-caustic deviations involve a higher degree of model 
degeneracies with more prominent finite-source and parallax ef- 
fects. In any case, a promising sensitivity to Earth-mass planets is 
only reached for lens-source impact parameters ito < 5 X 10 -4 , 
which occur at a rate of less than one per year. 

On the other hand, the non-negligible probability to detect 



planetary signals if the source passes in the vicinity of planetary 
caustics offers a fair chance of detecting a planet of Earth-mass by 
also making use of the large number of events that exhibit lower 
magnifications at a given time. Given these facts, it is not a surprise 
that the first sub-Neptune mass planet whose existence could be 
reported on t he basis of microlen sing observations, OGLE 2005- 
BLG-390Lb dBeaulieu et al]|2006l) , produced a 15 to 20 per cent 
signal at a magnification A ~ 1.3 about 10 days after an observed 
peak at magnification Ao ~ 3 (see Fig. [TJ rather than a deviation 
within a highly-magnified peak. 

While the mass of OGLE 2005 -BLG-390Lb is about 5 M®, 
uncertain to about a factor of two dDominikl l2006t) , a planet of 
1 M© in the same spot would still have produced a signal with 
an amplitude of ~ 3 per cent, lasting ~ 12 h rather than about 
twice that long. The actual sampling would have been insufficient 
for discovering such a planet in this configuration, but the situation 
would have been different had we decreased our sampling inter- 
val to 10-15 min on the suspicion of a first deviation. This case 
explicitly shows how an anomaly detector can help us in not miss- 
ing short-lasting small deviations (related to low-mass planets). By 
requiring an initial sampling that is just dense enough for an on- 
going anomaly being alerted before most of it has passed, it more- 
over allows to monitor a sufficient number of events for provid- 
ing a reasonable number of planet discoveries. The main gain of 
the anomaly detector will indeed be achieved for detecting planets 
from perturbations related to planetary caustics at lower and mod- 
erate magnification, whereas a high-cadence sampling can already 
be scheduled a-priori for (predictable) high magnifications without 
the need for any further alert. 

The ability of detecting an anomaly depends on how well ear- 
lier data constrain the model describing an ordinary light curve. For 
large model parameter uncertainties, it becomes hard to distinguish 
a real deviation from a necessary model revision due to a previous 
misestimate, for which x 2 adjustments are not a reliable indica- 
tor due to the intricate parameter space and poor knowledge about 
the measurement uncertainties. Therefore, the anomaly detection is 
more efficient afte r the peak of a microlensing has passed rather 
than prior to it (c.f. lUdalskill2003l) , where the ability is particularly 
vulnerable to data gaps. Thus, if the increased detection efficiency 
for low-mass planets that is achieved by means of the anomaly de- 
tector is a relevant goal for a monitoring strategy, it is sensible to 
give preference to events past peak over those pre peak for compa- 
rable magnifications. Although it is more difficult to decide whether 
a deviation from a previous model is real or due to a model mises- 
timate if constraints on its parameters are weaker, it is more likely 
that a suspected deviation occurs and is reported. This has the by- 
effect that more data will be collected in this case, which in turn 
strengthens the model parameter constraints. Despite the fact that 
the higher magnification around the peak allows for accurate data 
being taken with shorter exposure times, the weak constraints on 
the position of the peak make it rather difficult to detect an ongoing 
anomaly there, unless the peak region is monitored quite densely 
and no data gaps occur. 



4 CONCEPT FOR AN ANOMALY DETECTOR 

If reported data deviate from the expected light curve, this could 
either mean that there is a real effect, the deviation could be of sta- 
tistical nature, or the data could simply be erratic by any means. It 
is therefore impossible to arrive at an appropriate judgement about 
the presence of anomalies on the basis of a single deviating data 
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point. However, such a point should raise suspicion that an anomaly 
is indeed ongoing. Our anomaly detector, dubbed SIGNALMEN, 
profits from the fact that real-time photometry and robotic telescope 
operation allow immediate feedback. Rather than having to rely on 
a fixed sampling rate for a given event, we can request prompt fur- 
ther observations once the modelling of incoming data indicates a 
deviation from an ordinary light curve. 

Based on the collected data, the anomaly detector can arrive 
at one out of three possible conclusions and assign a corresponding 
status to the event: 

• there is no ongoing anomaly (ordinary) 

• there is an ongoing anomaly (anomaly) 

• not sure what is going on (check) 

While the last option, corresponding to a suspected, unconfirmed 
anomaly, does not look appealing at first sight, it actually marks 
the strength of the feedback concept. In this case, we urgently re- 
quest further observations on the same target, thereby providing the 
anomaly detector with further data on which it can base the decision 
in subsequent runs. In a 'recheck and repeat' strategy, data whose 
absolute model residual is among the largest trigger further obser- 
vations, and this process is repeated until a decision about whether 
there is an anomaly can be taken with the desired significance. 

The art of optimizing an anomaly detector is in finding the 
appropriate balance between not missing planetary anomalies and 
avoiding false alerts. The availability of immediate feedback opens 
the possibility of using a rather low initial trigger level on the first 
suspicion of an anomaly, which gives us a fair chance of detecting 
low-amplitude anomalies at an early stage. The early detection is a 
vital feature for being able to discover Earth-mass planets. In con- 
trast, we do not care that much about the detection of anomalies 
that have already been missed or are mostly over. A low initial trig- 
ger however means that we will need to spend a significant amount 
of time on collecting evidence against the presence of an anomaly 
if the point that triggered the 'check' observations does not consti- 
tute a real deviation. As pointed out in more detail in the following 
section, we aim at rechecking 5 per cent of the incoming data for 
anomalous behaviour, while about 4 to 5 further points are expected 
to be required for providing sufficient evidence against. This means 
that we spend about 20 per cent of our observing time on checking 
potential anomalies. By basing the criterion for a significant devia- 
tion on a comparison of the model residual of the tested data point 
with those of earlier data, we pay respect to the fact that the true 
scatter of data is not properly reflected by the size of the reported 
error bars and can be non-Gaussian. 

We also account for the fact that data collected with differ- 
ent telescopes may arrive in blocks rather than point-by-point and 
not necessarily in time sequence. Moreover, all data are subject to 
change, which not only means that reported (F^ctf;) might al- 
ter between two runs of the anomaly detector, but data at certain 
epochs might disappear, whereas additional data at other epochs 
prior to the most recent data point might be released. By not making 
any distinction between whether 'new' data are released in a block 
or arise from recent point-by-point observations, we also take care 
of the possibility that an anomaly is already apparent in the latest 
data update. 

Our robust fitting scheme is rather powerful in identifying out- 
liers and therefore gives us some protection against failures of the 
real-time photometry and weird results that might be the conse- 
quence. We have implemented a further test for distinguishing be- 
tween havoc photometry and ongoing anomalies which produces an 
alert urging to check the data reduction. However, there is no way 



getting around the point that the capabilities of the anomaly detec- 
tor will rise or fall with the quality of the real-time data analysis. 
In principle, one can also investigate correlations with observing 
conditions such as the reported seeing or sky brightness. However, 
such information may not be provided for all considered sites, so 
that we try to avoid relying on it as long as possible. 



5 ANOMALY DETECTOR EMBEDDING AND 
EXTERNAL INTERFACES 

The intelligent-agent architecture of the eSTAR project consti- 
tutes the harness inside which the SIGNALMEN anomaly detec- 
tor operates. Thereby, it provides autonomous decision-making by 
means of software, which allows to build systems that learn and 
adapt. The eSTAR system provides the feedback loop by feeding 
the SIGNALMEN anomaly detector with real-time data, which 
then replies with an expert opinion that allows the eSTAR sys- 
tem to solve the distributed-scheduling problem of how to dis- 
tribute follow-up requests over the network in order to maximize 
the chances of detecting and characterizing an extra-solar planet. 

The eSTAR project serves as a meta-network between exist- 
ing proprietary robotic te lescope networks b uilt upon a peer-to-peer 
agent based architecture (Wooldridge 2002), which cuts across tra- 
ditional notions that running such a network requires a "master 
scheduler". Instead, eSTAR can be viewed as a collaborative multi- 
agent system using a contract model. The crucial architectural dis- 
tinction of such a system is that both the software controlling the 
science programme and those embedded at the telescope acting as 
a high-level interface to the native telescope control software are 
equally seen as "agents". A negotiation takes place between these 
agents in which each of the telescopes bids to carry out the work, 
with the user's agent scheduling the work with the agent embed- 
ded at the telescope that promises to return the best result. This 
preserves the autonomy of individual telescope operators to imple- 
ment scheduling of observations at their facility as they see fit, and 
offers adaptability in the face of asynchronously arriving data. For 
instance, an agent working autonomously of the user can change, 
reschedule, or cancel queries, workflows or follow-up observations 
based on new information received. The eSTAR architecture repre- 
sents a "turn-key" system for autonomous observations of transient 
events, and therefore is ideal for microlensing follow-up. 

The agent s are also capable of responding in real tim e to 
external alerts jWilliams & Seamanll200a ; Iwhjte et alj|2006h , so- 
called Virtual Observatory Events (VOEventsfl While OGLE and 
MOA alerts are being translated into this format, the detection of 
an anomaly by SIGNALMEN will also be reported by means of a 
VOEvent. 

Besides the communication by means of software agents, the 
design of the SIGNALMEN anomaly detector also contains inter- 
faces for output to human observers and upload of data provided 
by any other observing campaign. Currently, data from PLANET, 
OGLE, MOA, and MicroFUN are fed in. Moreover, we will keep 
two separate mailing lists for notification on the decision in favour 
of an ongoing anomaly ('anomaly' status) and on the detection 
of deviant points ('check' status), which everyone is free to sub- 
scribe to. While dense follow-up by other teams is much encour- 
aged in this case, the 'check' status will be invoked frequently 
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(several times each night) and mainly serves to steer the inter- 
nal feedback with the robotic telescopes of the RoboNet network 
and in second instance with the other telescopes involved in the 
PLANET/RoboNet campaign. In addition to providing real-time 
notification of suspected or ongoing anomalies, we will publish up- 
to-the-minute plots showing collected data along with a model light 
curve, whose parameters have been determined by the SIGNAL- 
MEN anomaly detector. 

On the suspicion of an anomaly, a fast response with further 
observations is crucial for either confirming or rejecting this hy- 
pothesis. While robotic telescopes can react almost instantaneously, 
human observers need to be informed by e-mail or other means 
of communication, which adds some delay. Only if an observa- 
tory is staffed and the observer frequently monitors incoming e- 
mail, the feedback loop can be closed. This works reasonably well 
with the current PLANET network, where observers are present at 
the telescope on each night with suitable weather during the ob- 
serving season. However, telescopes that are only activated on a 
target-of-opportunity basis, such as several of those used by Micro- 
FUN, might miss the short-notice call. In any case, the success of 
the strategy is limited by the need to find out whether a suspected 
anomaly is present or not with the use of telescopes that have al- 
ready monitored the microlensing event of interest. The value of 
data from other sites is limited to providing early useful data if it 
turns out that an anomaly is ongoing, but these contain rather little 
information about whether the light curve deviates. 

While so far, we have implemented an algorithm that alerts us 
on suspected or ongoing anomalies, it neither gives us a recommen- 
dation of the best anomaly sampling interval, for which we simply 
assume an initial choice of 10 min, nor does it inform us when the 
anomaly is over and we can return to the standard follow-up sam- 
pling rate. Both of these issues currently need to be dealt with by 
human interaction through an internal webpage automatically list- 
ing events that are considered to deviate from ordinary light curves. 



6 THE ANOMALY DETECTOR ALGORITHM 

6.1 Basics, data statistics, and deviations 

The implementation of the SIGNALMEN anomaly detector de- 
scribed in the following is a first sketch, matching the primary re- 
quirements. It involves some basic statistical tests, building upon 
prior experience. More sophisticated tests can be designed and 
added, should it turn out that these yield significant improvements 
on the decision process. During the 2007 season, SIGNALMEN 
will log all incoming data, the anomaly indicators, current model 
parameters, and its decisions, which will provide a valuable basis 
for further tuning. Our algorithm involves several constants that can 
be adjusted. Their values can be changed by editing a configuration 
file rather than requiring alteration of the source code itself. In the 
following, we list our default setting in brackets. 

With the source and background fluxes, and , depend- 
ing on the data archive s, residuals need to be compared by means 
of the magnifications 



s(i)} 



«] 



(9) 



rather than the measured fluxes F», where the uncertainties of Ai 
are given by 

a Ai =a F J\F^\ (10) 



In general, the reported error bars a Fi are not a proper reflec- 
tion of the true scatter, which moreover frequently deviates from a 
Gaussian distribution. In particular, data provided by OGLE come 
with severely underestimated photometric uncertainties for J < 15, 
whereas these are about the right size for 15 ^ J ^ 18 and over- 
estimates for faint targets I ^ 18. One of the sources of this be- 
haviour is that the photometric reduction packages usually do not 
take into account further systematic uncertainties. We therefore cor- 
rect for this fact by adding a systematic error SYST_ERR (0.003) 
in quadrature to the uncertainty of the reported magnitude. More- 
over, rather than relying on a Fi , we assess the scatter by means of 
two statistics, namely the median scatter and the critical scatter 
^crlf ^ calculating the residuals 
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and with the critical real index 
v [s] = n M ( X _ DEV _PERC) + 1 
we determine 
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where j3 [s] =v [s] ~ [v [s] J . DEV.PERC (0.05) denotes the probability 
for an absolute deviation in excess of <5^J t , and n' s ' is the number 
of data points for archive s. With a deviation threshold DEV_SIG 
(2), we require for a significant deviation both 

\Aj - A(tj)\ > DEV_SIG<7.4 3 max {l, S ls(j)] } 



and 



\A, - A(t 3 



> a Aj O clit 



(15) 



(16) 



to hold. For Gaussian errors bars, both conditions can be made to 
coincide. In order to allow for a proper evaluation of the scatter, 
we require that at least MIN_DATA_TEST (6) data points and data 
from at least MINJJIGHTS (2) previous nights have been collected. 
Otherwise, the statistical scatter is likely to be underestimated and 
therefore false alerts are almost certain. 

With our robust-fitting algorithm that downweighs or even ig- 
nores outlier and the fact that we rely on the median scatter and 
trigger on the absolute residual exceeding that of a fixed percentage 
of data, we are well able to distinguish between low-quality data 
and real deviations. In particular, this allows us to achieve a low 
false alert rate. The requirement of obtaining significant data statis- 
tics before assessing deviations come at the price of some inability 
to identify deviations in fast-rising events with high-magnification 
peak. However, this does not significantly affect the planet detec- 
tion prospects, since a high-cadence sampling will be carried out 
for these events irrespective of suspected anomalies in the data. 



8 M. Dominik et al. 



6.2 Data sequence and modelling 

SIGNALMEN assumes that events do not exhibit anomalies at the 
time these are first announced by the OGLE or MOA microlensing 
surveys. For each data archive we keep track of the latest collected 
data point and restart our assessment for anomalies at the epoch 
t ne w corresponding to the earliest 'new' point among all archives. 
In order to assess the data point by point, we sort these in time se- 
quence and step through points k ^ n with tk ^ t nmr , where n is 
the index of the most recently observed data point. For each event, 
we store the time ranges for which anomalies were considered to 
be ongoing, and the parts of these ranges prior to t ncw are then ex- 
cluded from fits for models of an ordinary light curve. Moreover, on 
each run of SIGNALMEN on a specific event, we note the epoch 
t c 1? £ncw for which an ongoing anomaly was first suspected, and 
administrate a list of all points I with t\ ^ t c that were found to 
deviate, which form the current anomaly sequence. When consid- 
ering all data with t ^ tk, the deviation of a point with index j 
(t c ^ tj ^ tk) can be determined with respect to the following 
models that include all data with indices i that fulfill: 

• 'previous': ti < t^, exclude data within an anomaly time 
range prior to t Ilcw or in the current anomaly sequence 

• 'current': ti ^ tk, exclude data within an anomaly time range 
prior to t ncw or in the current anomaly sequence 

• 'all non-deviant': ti ^ tk, exclude data within an anomaly 
time range prior to t new or in the current anomaly sequence, but 
include i — j 

• 'all-data ':ti ^ tk, exclude data within an anomaly time range 
prior to t new 

If there is no current anomaly sequence, i.e. none of the points 
k ^ n has been found to deviate so far, the 'all-data ' and 'all non- 
deviant' models coincide with the 'current' model. Since model 
predictions can be expected to fail frequently, our initial assess- 
ment of a deviation is solely based on the 'current' model, which 
includes the latest considered point k. Should this point fail to de- 
viate significantly by means of the conditions given by Eqs. dl5t 
and dl6t . the 'current' model becomes the 'previous' model and k 
is increased. Otherwise, t c = tk and data point k is added to the 
current anomaly sequence. While the 'previous' model is retained, 
it also becomes the 'all non-deviant' model, whereas the 'current 
model' also becomes the 'all-data' model. For increased k, further 
tests will be performed for data j (t c ^ tj ^ tk)- 



6.3 Anomalies: accept or reject? 

If a current anomaly sequence has been found, SIGNALMEN will 
try to figure out whether further data points provide evidence in 
favour of an ongoing anomaly or against it, leading to finishing up 
with 'anomaly' or 'ordinary' status. If the current data do not allow 
to arrive at either of these conclusions, the 'check' status is invoked. 
In this case, the markers for the latest data points for each of the 
archives are set so that the current anomaly sequence is reassessed 
on the next run of SIGNALMEN. This avoids the necessity to 
store further information about this sequence and also easily allows 
for a potential revision of these critical data in the meantime. 

Data taken after t c that are found not to deviate significantly 
from the 'current' model can provide evidence against the presence 
of an ongoing anomaly. However, simply counting the number of 
non-deviant points is not a viable option since these might have 
larger uncertainties than the deviant points. This happens in partic- 
ular if later data originate from different sites, while even for the 



same site it cannot be guaranteed that the same photometric uncer- 
tainty can be retained. Since data with large scatter and therefore 
no indication of an anomaly must not be used as evidence against, 
it is unavoidable that the photometric uncertainties are taken into 
account. Moreover, we also need some characteristic for the ampli- 
tude of the anomaly which we would like to decide about whether 
it is real or not. Let us consider the fractional deviation 

£l ~ A{U) ' U/) 
and for a deviant point I define ei as the anomaly level. With a Ej — 
(o'Aj max 1 1, 5^"^\)/A(tj), we then obtain the weighted aver- 
age over all non-deviating points j after the deviant point (i.e. tj > 
ti) 
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and its standard deviation 
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The anomaly is then called off if 
e < ei/2 - REJECT.SIG a- 
£ > £i/2 + REJECT.SIG a T 



(for ei > 0) 
(for ei < 0) 



(18) 



(19) 



(20) 



with a default setting REJECT.SIG = 1 and the additional require- 
ment that at least MINPTS.REJECT (4) non-deviant points have been 
collected. For Gaussian data with constant error bars, we find the 
anomaly call-off typically not requiring more than 5 measurements. 
However, this can take significantly longer if only data with large 
effective error bars (corrected for actual scatter) can be acquired. 

If the data point k has been found not to deviate, we also re- 
assess the current anomaly sequence with respect to the 'all non- 
deviant' model. If an anomaly cannot be confirmed or discarded, 
just testing points in sequence against the 'current' model can ei- 
ther lead to missed anomalies or false alerts if the model is not 
well-constrained. We therefore determine the residuals with respect 
to a model that includes all points found not deviating (and their 
scatter). This also allows taking into account an increased scatter 
present in more recent data. Previously deviant data that do not 
fulfill the new criterion are removed from the current anomaly se- 
quence, which might lead to a revision of t c and leave SIGNAL- 
MEN with an empty current anomaly sequence. In the latter case, 
SIGNALMEN will continue as if no deviant points were found in 
the current run. We also require that all data points in the current 
anomaly sequence deviate to the same side. Therefore, it will be 
shortened if necessary to meet this condition. 

Similarly, if the most recently considered data point k is found 
to deviate to the opposite site as the previous data, a new current 
anomaly sequence is started at t c = tk and the previous sequence 
is abandoned. 

A stronger hint for an anomaly being ongoing is obtained if 
the data point k deviates to the same side as the previous points in 
the current anomaly sequence. Once the current anomaly sequence 
contains at least two data points, we start testing the collected data 
against an 'all-data' model, which also contains the points in the 
current anomaly sequence. With the earlier tests we avoided that 
the model of an ordinary event is driven towards points that devi- 
ate from it, which allows us to call off an anomaly if further points 
follow an ordinary light curve without getting confused by outliers. 
However, we also need to take care of the fact that more weight than 
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just that of a single point might be needed to correct for a bad ear- 
lier estimate of model parameters. As a compromise, we adopt less 
strict criteria, namely that the residuals of the last MINPTS.ANOMALY 
(5) points are all of the same sign and at least MINPTS_ALL_ANOM 
(3) points deviate significantly. If earlier data in the current anomaly 
sequence cannot match these criteria, the sequence is shortened and 
t c is revised. 

A further test is based on the overlap between the points in 
the current anomaly sequence and non-deviant points falling in be- 
tween. With the 'all-data ' model, we determine 

SAi =Ai- A(ti) . (21) 

If for a non-deviant point j following a deviant point / for which 
SAi > 0, one finds 

SAj + 2o A] max{l,<5 [s(j)1 } < SA t - 2a Al max{l, <5 [s(01 } (22) 

or the equivalent relation to hold for the subsequent deviant point, 
the non-deviant point is considered to contradict point / deviating, 
which is therefore removed from the current anomaly sequence. For 
8 A t < 00 the corresponding condition reads 

SAj ~ 2a A] max{l,5 [s(j)1 } > SA l +2a Al max{l, S ls(l)] } .(23) 

Finally, we realize that the photometric reduction might fail 
occasionally and produce weird results. A common characteris- 
tic that can be distinguished from real anomalous behaviour are 
sudden changes between a rise and fall. We therefore determine 
the pattern of significant increase or decrease of the magnification 
amongst the data in the current anomaly sequence. Should there be 
more than one change in direction, SIGNALMEN abstains from 
the claim that an anomaly is ongoing. This 'zig-zag test' is only 
used as the final criterion once all other conditions for an ongo- 
ing anomaly are fulfilled. For two deviant points I and m > I, a 
significant increase is characterized by 

SA m -2a Am max{l,<5 [s(m)1 } > SA t +2a Al max{l, <5 [s(01 } , (24) 

whereas a significant decrease arises by exchanging I and m. If 
there is no significant change between neighbouring points, a sig- 
nificant increase is assessed with respect to the lowest of these 
points while a significant decrease refers to the highest of these. 

To summarize, SIGNALMEN concludes that there is an on- 
going anomaly if all of the following criteria are satisfied 

• the anomaly is not called off by means of a series of at least 
MINPTS_REJECT (4) non-deviant points with a weighted-average 
fractional deviation significantly (measured by REJECT_SIG (1.0)) 
closer to zero than half of the fractional deviation of the previous 
deviant point 

• the most recent deviant points form a sequence of at least 
MINPTS_ANOMALY (5) points that were found to deviate to the same 
side from the 'current' model and the 'all non-deviant' model 

• the residuals with respect to the 'all-data ' model of at least the 
last MINPTS_ANOMALY (5) points in the current anomaly are all of 
the same sign 

• at least MINPTS_ALL_ANOM (3) points in the current anomaly 
sequence deviate from the 'all-data ' model 

• there are no non-deviant data points in between those in the 
current anomaly sequence that significantly fail to overlap with 
them 

• data in the current anomaly sequence do not change more than 
once between a significant increase and decrease 

8 Obviously, there is no SAi = case. 



If these criteria are fulfilled for k — n, i.e. at the end of 
the collected data, SIGNALMEN activates the 'anomaly' mode. 
Should these be fulfilled earlier (k < n) only, SIGNALMEN 
finishes with 'ordinary' status, but a file notifying about a missed 
anomaly is written. If just the zig-zag test fails, SIGNALMEN no- 
tifies about problems with the photometric reduction and suspends 
evaluation of data archives for which recent data showed more than 
one change of direction in the suspected anomaly sequence. Such 
a case needs human intervention and should be dealt with at high 
priority. 



7 PROSPECTS WITH THE ANOMALY DETECTOR 

In order to demonstrate what can be achieved with the anomaly de- 
tector, let us use the event OGLE 2005-BLG-390, which already 
allowed us to detect a planet of 5 M® (with a factor two uncer- 
tainty), as an illustrative example and starting point of the dis- 
cussion. Fig. [2] shows the model light curve for the correspond- 
ing configuration again, where the planet OGLE 2005-BLG-390Lb 
has been replaced by a 1 M® version in the same spot, but rather 
than the collected data, we now show simulated data related to the 
three robotic 2m telescopes that currently comprise the RoboNet- 
1.0 network: the Liverpool telescope (LT), the Faulkes telescope 
North (FTN), and the Faulkes telescope South (FTS). According 
to the target observability from the different sites at the event co- 
ordinates (RA = 17 h 54™19P19, Dec = -30?22f38'.'3 (12000)), 
requiring that the target is at least 20° above the horizon, synthetic 
data have been generated where the average sampling interval is 
At = (2 h)/y/A and the photometric accuracy is 1 per cent at 
baseline and smaller as the event brightens, following photon noise 
statistics, where Gaussian errors have been assumed. A systematic 
error of 0.5 per cent has been added in quadrature. For the time 
of the next observation, a Gaussian fluctuation of 20 per cent of 
the sampling interval has been adopted, and while its photometric 
uncertainty itself fluctuates by 12.5 per cent. Moreover, a drop-out 
probability of 5 per cent on each data point has been assumed. 

While it would have been rather easy to detect and character- 
ize a planet like OGLE 2005-BLG-390Lb, the standard sampling 
would have given rather little evidence for an Earth-mass planet in 
the same spot, and a characterization would not have been possi- 
ble. Arrows in Fig. [2] indicate data points that were found to devi- 
ate by the anomaly detector, given the best-fitting model that could 
have been obtained at that time, based on all previous data. Fur- 
ther data after the first four trigger points would have indicated that 
there is no ongoing anomaly, but the deviation is rather of statisti- 
cal nature. In contrast, subsequent data points after the first trigger 
point that falls onto the real anomaly would have confirmed the de- 
viation and finally led to the activation of high-cadence anomaly 
monitoring. This example however also shows a critical weakness 
of the current RoboNet-1.0 network, namely its lack of round-the- 
clock coverage. In particular, one sees that the southern telescope 
offers a much longer time window for our purpose than either of the 
northern telescopes. Just after the opportunity of taking a further 
point after activation of the 'check' status, the target could not have 
been followed anymore with the FTN and it would have been nec- 
essary to wait for the LT for acquiring subsequent measurements. 
This demonstrates that provision of a fast response also implies the 
availability of a telescope at the requested time. Further telescopes 
available in South Africa and Chile (see Fig.[2jl would have allowed 
a coverage of the anomaly sufficient to detect an Earth-mass planet. 

The detection of OGLE 2005-BLG-390Lb was eased by the 
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OGLE 2005-BLG-390 -> Earth-mass planet OGLE 2005-BLG-390 -> MS star, m = 1M, 
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Figure 2. Detection of a hypothetical Earth-mass planet in microlensing 
event OGLE 2005-BLG-390 located in the same spot as OGLE 2005-BLG- 
390Lb with the RoboNet telescope network. Simulated data for the differ- 
ent telescopes are shown in different colours along with theoretical light 
curves, where a thin solid line corresponds to the actually derived model 
for that event (with OGLE 2005-BLG-390Lb), a bold solid line to a model 
with an Earth-mass planet, and a bold dotted line to a model without planet. 
Aitows mark data points that have been found to deviate from the best- 
fitting model available at that epoch. While the top panel shows only data 
collected with the standard sampling, the middle and bottom panel include 
further data with high-cadence (10 min) sampling after having triggered on 
the anomaly. For the top and middle panel, only the current RoboNet- 1.0 
telescopes have been considered, whereas for the bottom panel, the avail- 
ability of two additional similar telescopes in Chile and South Africa has 
been assumed. 



source being a giant star with a radius i?* ~ 9.6 Rq, which not 
only allowed to obtain an accurate photometry with rather short 
exposures but also increased the probability of detecting a signif- 
icant deviation. While the large angular size led to a reduction of 
the amplitude, which did not matter because it remained at the 15 
to 20 per cent level, a larger range of orientations or impact param- 
eters than for a smaller source would have created on observable 
signal. Moreover, the duration of the planetary deviation is domi- 



7 1 


1 1 1 1 | 1 



















10.5 11 



• LT 

FTN • Chile 

•FTS S Africa simulation 

_, I I I I I I I I L_ 



5 10 

days since 31.0 July 2005 UT 

Figure 3. Detection of a hypothetical Earth-mass planet in an event resem- 
bling OGLE 2005-BLG-390 at the same separation as OGLE 2005-BLG- 
390Lb, but for a main-sequence source star with i?* ~ 1 . 2 Rq . In fact, the 
source size has been assumed to be 8 times smaller than for the giant ob- 
served in OGLE 2005-BLG-390. Since with the smal ler source, t he planet 
is missed for the oris inal angle a = 157.9° teeaulieu et alj2006h between 
the line from the planet to its host star and the source trajectory, where the 
lens centre of mass is to the right, a slightly different angle a = 158.2° 
has been adopted, for which a 5 per cent deviation results. Otherwise, this 
figure is analogous to Fig. [2] 



nated by the source star moving by its angular size, giving a rather 
comfortable timespan, which would still have been ~ 12 h for an 
Earth-mass planet. While a main-sequence star could have provided 
a signal with larger amplitude, the probability to observe it would 
have been smaller and it would not have lasted for that long. 

If one replaces the source star of OGLE 2005-BLG-390 with 
an 8 times smaller version, an Earth-mass planet in the same spot 
as OGLE 2005-BLG-390Lb would become undetectable since the 
smaller source would not enter a region for which significant de- 
viations result. For the configuration shown in Fig. [3] the angle of 
the source trajectory relative to the planet-star axis has therefore 
been slightly adjusted, resulting in a 5 per cent deviation. Achiev- 
ing good photometry on the fainter target is more difficult and re- 
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Figure 4. Coverage of a microlensing event resembling OGLE 2005-BLG- 
390 with a main-sequence source star (R* ~ 1.2 Rq) instead of the 8 
times larger giant, and a planet of 0.1 Mj at an angular spearation of d = 
1.25 times the angular Einstein radius 6>e- The angle between planet-to-star 
axis and source trajectory, with the lens star to the right, has been chosen 
as a = 126.6°, so that the source crosses the planetary caustic. As for 
Fig. [2] the trigger points of the anomaly detector are indicated by arrows, 
and the middle and bottom panels include the high-cadence follow-up after 
anomaly suspicion or detection. While the top and middle panels only show 
data corresponding to the current RoboNet- 1 .0 network, the availability of 
further similar telescopes in Chile and South Africa has been assumed. 

quires longer exposure times. Nevertheless, PLANET has demon- 
strated a photometric accuracy of even less than 0.5 per cent on 
a main-sequence star is possible provided that it is fairly isolated 
rather than in a crowded area. While for the previously discussed 
case involving a giant source star, signal amplitudes significantly 
exceeding 3 per cent cannot result, the shown 5 per cent devia- 
tion is not even near the limit for main-sequence stars, for which 
very strong signatures become possible should the source happen 
to cross a caustic. One also sees that the duration of the planetary 
deviation has not decreased by a factor of 8 as the source size did. 
Contrary to the giant source star case, the signal duration is now 



roughly given by the time in which a point source passes the re- 
gion of angular positions that lead to significant deviations, and 
remains ~ 8 h for the prominent peak. The angular size of the 
source star itself is reflected in the small peak within the timespan 
over which the brightness in presence of the planet is smaller than 
without. As before with a giant source star, the proper characteri- 
zation of the planetary anomaly is not possible with the standard 
2 h sampling, while high-cadence sampling after having suspected 
or detected an anomaly will solve the problem, provided that tele- 
scopes are available to observe the target. Interestingly, in a very 
early stage of the anomaly, one of the data points appears to be 
higher just by chance. Further data taken at the sampling interval 
of 10 min however do not confirm a significant deviation, so that 
only after the next data point taken with the standard sampling rate 
the high-cadence anomaly monitoring remains active. 

After having found that the discovery of Earth-mass plan- 
ets does not constitute the limit of what can be achieved with 
microlensing survey/follow-up campaigns equipped with an auto- 
mated anomaly detector, let us look into how far one can go. In 
fact, the rather large separation d ~ 1.6 of OGLE 2005-BLG- 
390Lb from its host star did not offer a very fortunate configu- 
ration. Let us therefore also consider d ~ 1.25 and see how the 
signal amplitude and duration are affected. As Fig. [4] shows, even 
for a planet with mass m — 0.1 M® located at 1.25 8e from its 
host in an OGLE 2005-BLG-390-like event, a signal of 10 per cent 
lasting about 2.5 hrs can result on a main-sequence source star. The 
early detection of such a short signal with a standard survey sam- 
pling interval of 1-2 hrs and the anomaly monitoring each become 
challenging. In fact, the rather short time gap of about 40 min be- 
tween the FTN in Australia and a telescope in South Africa is suf- 
ficient for missing the crucial falling part of the planetary anomaly 
in our simulation. This also demonstrates the extraordinary value 
of a telescope at the western edge of Australia and/or in southern 
central Asia. Nevertheless, the discovery of planets with masses of 
even 0.1 or below by ground-based microlensing campaigns 
remains a matter of probability rather than possibility. 



8 SUMMARY AND CONCLUSIONS 

For probing models of planet formation and evolution and thereby 
taking an important step towards our understanding of the origin 
of living organisms as we know them, microlensing will remain 
a competitive and complementary technique with respect to other 
methods for the foreseeable future with its unique sensitivity to 
low-mass planets in wider orbits. It is unlikely (although not im- 
possible) that microlensing will provide a timely discovery of a 
planet on which conditions similar to those on Earth can exist that 
are known to support the formation of life forms, with less than 
3 per cent of planets in a ny mass range ex pected to orbit suitable 
host stars at suitable radii jParketall 2006). While both transit and 
radial-velocity surveys approach the required orbital range for such 
habitable planets from closer orbits, essentially all planets that can 
be expected to be detected by microlensing reside in wider orbits. 
The discovery of the first extra-solar planets already demonstrated 
impressively how little one can infer about the origin of the Solar 
system if the study remains restricted to itself. Similarly, one should 
not expect that a study just of habitable planets will allow to arrive 
at a well-understood picture of their formation. Instead, a reliable 
test of theories should involve data spanning over a wider and sur- 
rounding region. Moreover, microlensing allows to obtain a planet 
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sample not only around stars in the Galactic disk but also around 
those the Galactic bulge, thereby probing two distinct populations. 

We have shown that our SIGNALMEN anomaly detector, 
which went into operation for the 2007 microlensing observing 
season, allows to adopt a 3-step strategy of survey, follow-up, and 
anomaly monitoring. The basis of this strategy is formed by the 
microlensing events provided in real-time by the OGLE and MOA 
surveys out of which a sufficient number are then monitored with 
a network of 2m telescopes so that deviations due to Earth-mass 
planets are unlikely to be missed and a reasonable number of low- 
mass planets is expected be detected over the next few years. In 
particular, it is only required that the follow-up observations pro- 
vide an early enough trigger of an ongoing anomaly, whereas a 
proper characterization need not be ensured by these, because this 
will be achieved by the high-cadence anomaly monitoring after an 
anomaly has been suspected or detected. In fact, the use of an 
anomaly detector becomes quite efficient if many points are re- 
quired for proper characterization of a signal rather than being able 
to claim a detection from a single deviant point. The expected de- 
tections will provide a powerful test of models of planet formation 
and evolution around K- and M-dwarfs. While planets of Earth 
mass appear to have some particular appeal, they do not provide 
the hard limit for ground-based microlensing searches. As shown 
by one of the examples discussed in Sect. [7] the anomaly detector 
allows us to go even below 0.1 Mq, although such detections are 
not very likely to occur. However, these are reasonable goals for 
a network o f ground-based w ide-field telescopes or a space-based 
telescope ( Ben nett et ai]|2004l) . 

By only checking for significant deviations from an ordinary 
microlensing light curve, our current anomaly detector is blind to 
the nature of origin of the observed deviation. While it is ensured 
that planetary deviations due to planets of Earth mass and even be- 
low can be detected, more than 90 per cent of all deviations are due 
to other causes, such as finite-source effects, Earth-sun parallax, 
or stellar lens or source binaries. In order to distinguish these from 
planetary deviations and at the same time to obtain appropriate esti- 
mates for the urgency and frequency of second-level follow-up ob- 
servations, a full real-time modelling taking into account all these 
effects would be required. Optimally, the prioritization of events 
would follow the expected constraints on the planet characteristics 
rather than just maximizing their detectability while ignoring the 
chances of properly characterizing a potential planet. We plan to 
implement such a system in the near future. 
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Table Al. Forbidden regions for the planet separation as function of the 
event impact parameter 
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d± are the critical separations for which the centre of the planetary caus- 
tic(s) falls inside a circle of radius Mo around the host star. For sufficiently 
small mass ratios, this will prevent planets with separation parameters in the 
range (d— , d+ ) from being detected in events with impact parameter mo ■ 

APPENDIX A: THE INABILITY TO DETECT PLANETS 
AT THE RESONANT SEPARATION IN SOME EVENTS 

In general, angular positions for the source star relative to the lens 
composed of a star and its planet, for which a significant devia- 
tion in the observed light curve as compared to a lens star without 
planet results, form regions around the caustics of the lens sys- 
tem. For point sources, these regions grow as the angular separa- 
tion parameter d, where d0E gives the angular separation of the 
planet from its ho s t star, approaches unity. This indeed implies, as 
iMao & Paczvfiskil dl99lh pointed out, that planets around d ~ 1 
are most easily detected among all events that occur. However, this 
does not imply that this also holds for each of the observed events. 

Finite-source effects may cause an increase of the detection 
efficiency by means of the larger source catching more easily a 
region of significant deviation without bringing the signal ampli- 
tude below the detection threshold, or these can lead to a reduction, 
in particular if the finite source sub tends regions correspon ding to 
deviations of opposite signs. In fact lGaudi & Sacked l l2000h found 
that a reduction of the detection efficiency most prominently affects 
separations around the resonant d ~ 1, where the wide-separation 
side suffers more than the close-separation side. 

In any case, there is already another effect that prohibits the 
detection of less massive planets around d ~ 1. Given that the cen- 
tral caustic is found at the position of the lens star and the centre of 
the planetary caustics is separated by \d— l/d\ from it, the regions 
for which a significant deviation results might fall inside a circle 
whose radius is given by the impact parameter uo of a given event. 
In this case, these cannot be traced by the source trajectory. For 
small planet-to-star mass ratios q, the critical separations approach 

d ± = \ (V M o + 4 ± M o) (Al) 

~ 1 ± H (wo « 1) , (A2) 

while for larger q, an increasingly larger size of the caustics and 
the associated regions of significant deviations allows to enter the 
range d G (d_ , d+ ) further and further. The fact that for larger lens- 
source separations, less massive planets cannot be detected in in- 
creasingly broader regions around the angular Einstein radius #e is 
apparent in the examples showing the detection efficiency as func- 
tion of the separation parameter d for a few choi ces of the planet- 
to-star mass ratio in Figs. 5 and 8 provided by iGaudi & Sacked 
(2000), but unfortunately not discussed there. 

Table I All shows the critical separations for selected values of 
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uo- In particular, (low-mass) planets within the lensing zone can 
only be detected in events with uo ^ 1 (corresponding to Ao ^ 
1.34). Given these findings and the fact that a planet actually did 
reveal its existence in the event OGLE 2005-BLG-390, it is less 
surprising that with uo — 0.359, its angular separation from its host 
star is the rather large d — 1.61 angular Einstein radii, whereas a 
detection at e.g. d = 1.1 would have been impossible. 
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